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Abstract 

Charge density waves are thought to be common in two-dimensional electron systems in quantizing magnetic fields. Such 
phases are formed by the quasiparticles of the topmost occupied Landau level when it is partially filled. One class of charge 
density wave phases can be described as electron solids. In weak magnetic fields (at high Landau levels) solids with many 
particles per unit cell - bubble phases - predominate. In strong magnetic fields (at the lowest Landau level) only crystals 
with one particle per unit cell - Wigner crystals - can form. Experimental identification of these phases is facilitated by the 
fact that even a weak disorder influences their dc and ac magnetotransport in a very specific way. In the ac domain, a range 
of frequencies appears where the electromagnetic response is dominated by magnetophonon collective modes. The effect of 
disorder is to localize the collective modes and to create an inhomogeneously broadened absorption line, the pinning mode. In 
recent microwave experiments pinning modes have been discovered both at the lowest and at high Landau levels. We present 
the theory of the pinning mode for a classical two-dimensional electron crystal collectively pinned by weak impurities. We 
show that long-range Coulomb interaction causes a dramatic line narrowing, in qualitative agreement with the experiments. 
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1. Introduction 

Two-dimensional electron gas (2DEG) in quantiz- 
ing magnetic fields is known to be a system with a 
rich phase structure. Recently, much attention has 
been devoted to charge-density wave phases, which 
include stripe phase, bubble phase, and a Wigner 
crystal [1,2,3]. An important new information on the 
dynamics of such phases was provided by a group of 
microwave experiments [4,5,6,7,8]. They demonstrated 
that the low-frequency ac response of magnetic-field- 
induced charge-density waves is strongly affected by 
disorder, in ways not anticipated in prior theoretical 
work on the subject. Motivated by these intriguing 
results, Huse and the present author [9] reconsidered 



the problem of collective dynamics of a Wigner crys- 
tal pinned by quenched disorder. Here I give a brief 
account of this work and its more recent extensions. 
Much of the foregoing discussion also applies to the 
bubble phase. The rest of this section is devoted to a 
brief introduction to the charge-density wave phases, 
conventional expectations regarding the response of 
a pinned system, and how they compare with the 
observed behavior. Sections 2 and 3 outline our theo- 
retical model and its analysis. Finally, Sec. 4 contains 
comparison of our theory with the experiment and 
concluding remarks. 

The bubble phase [3] can be viewed as a crystal made 
of multi-electron droplets (bubbles). The number of 
particles per bubble M depends on the Landau level 
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filling fraction v and changes in discrete steps at cer- 
tain critical values of v. Here M counts only the quasi- 
particles of the topmost (JVth) Landau level, which is 
partially filled. The other Landau levels are completely 
filled and in a first approximation, provide an inert 
background of constant uniform density. Bubble phases 
with M > 1 appear when three or more Landau levels 
are populated, v > 4 (N > 2). The Wigner crystal is 
the M = 1 bubble phase (one particle per unit cell) . It 
forms at any Landau level N provided v is close enough 
to an integer. We refer to the bubble and the Wigner 
crystal phases jointly as electron crystals. At zero tem- 
perature, in the absence of disorder and commensu- 
rability effects such crystals have a perfect triangular 
lattice and a long-range order. They possess collective 
modes referred to as magnetophonons. Wigner crystal 
has a single magnetophonon branch [10] , bubble phases 
may have several [11]; however, in both cases only one 
mode is gapless. The following discussion is devoted 
exclusively to the dynamics of the gapless mode. Qua- 
siclassically, a magnetophonon excitation can be un- 
derstood as a coherent precession of M-electron bub- 
bles around their equilibrium positions, with the charge 
distribution (form-factor) of individual bubbles being 
fixed. In this picture the dynamical degrees of freedom 
are the centers of mass of the bubbles. 

The effects of disorder on gapless magnetophonon 
excitations belong to a broad class of phenomena 
known as "pinning. The pinning affects a variety of 
physical systems. It has been extensively studied in the 
context of conventional charge-density waves [12] and 
vortex lattices in superconductors [13]. Two pinning 
mechanisms are distinguished, the individual and the 
collective ones. In this paper we focus on the latter. 
The collective pinning is the regime where the random 
potential due to impurities and other defects is too 
weak to significantly deform the crystalline lattice, so 
that the crystal is well ordered at small length scales. 
However, the cumulative effect of the disorder even- 
tually dominates the crystal elasticity at a pinning 
length, which is much larger than the lattice constant 
a. The long-range order of the crystal is thereby de- 
stroyed. Similarly, the low-frequency collective modes 
are unaffected by disorder at the scale of a, but get 
localized at some larger scale. These localized modes 
give rise to a maximum at a disorder-dependent fre- 
quency uj p in the real part of diagonal conductivity 
Keer(aj). Such maxima (pinning modes) have been 



studied primarily in conventional charge-density wave 
materials [12,14]. However, a few observations of simi- 
lar absorption lines due to 2DEG subject to a strong 
magnetic field have also been reported [15]. The afore- 
mentioned latest experiments on high-mobility 2DEG 
provided much more extensive and higher resolution 
data in the strong-field regime [4,5,6,7] and also de- 
tected pinning resonances in weak fields [8], which 
were interpreted as signatures of bubble phases and 
Wigner crystals at higher Landau levels. The quanti- 
ties that can be extracted and analyzed most reliably 
are (i) the position uj p of the peak in Ke a, (ii) its width 
Au p , and (iii) the area F under the curve (see Fig. 1). 
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Fig. 1. The sketch of the absorption line due to pinning. 

(i) The dependence of ui p on magnetic field and dis- 
order is often discussed in the framework of a popular 
physical picture [16,17] illustrated by Fig. 2. Each lo- 
calized magnetophonon mode is visualized as preces- 
sion (or drift) of a crystallite with linear dimension of 
the order of the pinning length R c in a potential well 
created by disorder. Let ui c = eB/m e c be the cyclotron 
frequency in the external magnetic field B. For a po- 
tential well with curvature ~ MqUq, where Mo is the 
total mass of the patch, an elementary calculation gives 
the frequency of the drift motion uj p ~ luq/uj c in the 
experimentally relevant case uj c ^> uj p . If u>a depends 
only on the amount of disorder (e.g., impurity concen- 
tration) but not on B, then ui p oc 1/B. 

(ii) As far as Atu p is concerned, in a random system of 
independent oscillators we expect a broad distribution 
of u)q [17], which implies a broad spectrum of their 
eigenfrequencies, i.e., Au p ~ u> p . 

(iii) Finally, for the integrated oscillator strength F, 
this model predicts [16], F = (it /2)(e 2 n e /m e )(u} p /uj c ), 
where n e and m e are the electron concentration and 
effective mass, respectively. 

Despite a certain appeal of these arguments, their 
predictions are not supported by experiments. First, 
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Fig. 2. Naive cartoon of a localized magnetophonon mode. 
The ball represents a single domain of the pinned crystal, 
which has mass Mo- The paraboloid symbolizes the pin- 
ning potential in which this domain resides. The dashed 
line shows the drift trajectory in the two-dimensional plane 
r = (x, y). 

the pinning frequency shows a complicated nonmono- 
tonic dependence on B. Second, the absorption line is 
very narrow. Quality factors Q — ui p /Acu p as large as 
ten were reported in the published literature on the 
high-field Wigner crystal [4,6,7] and lower but still re- 
spectable quality factors Q ~ 3 for the bubbles. Only 
the integrated oscillator strength F is in a rough agree- 
ment [8] with the simple theory. 

It turns out that the anomalous B-dependence can 
in principle be explained by the interplay between 
the correlation length of the pinning potential and 
the B-dependent form-factor of the electron wave- 
functions [18,19,9]. On the other hand, to explain the 
small linewidth one has to develop a theory that goes 
beyond the cartoon depicted in Fig. 2. Such a the- 
ory is the subject of this paper. We will show that 
the line narrowing is due to the long-range Coulomb 
interaction in the electron crystal. 



2. Definition of the model and the pinning 
frequency 

We treat Wigner and bubble phases on the same 
footing, as crystals of classical particles of mass Mm e , 
charge Me, and average concentration, n e /M. They in- 
teract with each other via potential V(r) and are sub- 
ject to a weak random potential U (r) . We assume that 
V and U incorporate the appropriate form-factors, cal- 
culated using the quantum wavefunctions of the bub- 
bles [1,3]. (This is the only place where quantum me- 
chanics enters our model.) We describe the deviation of 
the crystal lattice from the ideal periodicity in terms of 
the elastic displacement field u e ; = u' ' + u, which we 



split into the static ground-state distortion u' ' (r) and 
the magnetophonon fluctuations u(r, £). We restrict 
ourselves to the harmonic approximation, where mag- 
netophonon modes are obtained by diagonalizing the 
dynamical matrix [D] _1 = [D°]~ 1 <5 qii q 2 +<5S(qi — q2), 
where D° = D°(q, lj) is a dynamical matrix of a uni- 
formly pinned crystal. In the convenient basis of trans- 
verse (T) and longitudinal (L) components, u(q) = 
qiiL(q) + [z x q]uT(q), where q = q/|q|, D° has the 
form [16,10] 

rryji-i _ M 2 + So-pu 2 -ipuju c 

' ' ~ 2 2 ' 

ipcuuic Xq + So — P<j-> 

where p — m e n e is the average mass density, [i ~ 
M 2 e 2 /fta 3 is the shear elastic modulus, A = X(q) = 
2-Ke 2 n 2 / ' nq is the effective bulk modulus, and k is the 
dielectic constant of the medium. The q-dispersion of A 
is due to long-range Coulomb interaction. At small q, 
we have X(q) ^> p, i.e., compressions cost much more 
energy than shear deformations. Parameter So is the 
average value of the diagonal components of the cur- 
vature tensor S = VVC(r) evaluated at the ground- 
state particle positions. From standard collective pin- 
ning arguments, we expect So ~ fJ-/Rc = P^Oi where 
R c is the pinning length related to the amplitude of 
the random potential U and the shear modulus p [9]. 
The quantity 8S = S — Sol describes fluctuations of 
the curvature, where I is the identity matrix and tilde 
denotes the Fourier transform. Our goal is to calculate 
the ac conductivity, cr(u>) = —ie 2 n 2 uj{D(uj+iO)}, where 
(...) denotes disorder averaging. The conductivity can 
be expressed in terms of the self-energy II(q, u) = 
Sol + ((D)- 1 - [D ]" 1 ) , as follows: 

»e,H=e 2 ^9m--fHM_i . (2) 

[11(0, e) - e] 2 - ee c 

Here we introduced convenient "energy" variables e = 
pui 2 and e c = pui 2 . From Eq. (2) one can see that in 
strong magnetic fields the conductivity as a function of 
lo has a maximum of width Au) p = —Sim 11(0, pulp) / pui c 
centered at uj p — KeII/pa> c . Up to logarithmic fac- 
tors [9], 3?eII ~ So, and so ui p = ui 2 /ui c . Thus, the 
position of the conductivity peak is simply related to 
the pinning length R c , in basic agreement with the ar- 
guments given in Sec. 1 (and as mentioned above, R c 
can be straighforwardly expressed in terms of the mi- 
croscopic parameters of the model) . The lineshape and 
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linewidth require much more elaborate analysis given 
in the next section. 



3. Lineshape and linewidth of the pinning mode 

To find II and a we have to resort to certain approxi- 
mations, which we motivate by physical arguments. As 
a starting point, consider a simpler model of a pinned 
one-dimensional (ID) electron crystal with short-range 
(screened) interaction V in the absence of any mag- 
netic field. It has been known from early work [20,21] 
that 5ftecr(u;) has a broad maximum at the character- 
istic pinning frequency LUp~ D — v/R c , where R c is the 
pinning length and v is the sound velocity. However, 
there has been a disagreement about the behavior of a 
at small u>. Recently, the present author [22] and later 
also Gurarie and Chalker [23] advanced analytical and 
numerical arguments that $tea(ui) oc lu 4 for u> <C uJp D ■ 
The physical picture is as follows [22]. In ID, elastic 
displacement u is a scalar and there is a formal analogy 
between the phonon problem and ID localization in a 
random potential S(x) — U"(x). Indeed, the phonon 
mode with "energy" e; must satisfy the Schrodinger 
equation —pv 2 u" + S(x)ui(x) = e»M. Borrowing stan- 
dard arguments from the ID localization literature, one 
can show that phonon eigenmodes with lo < io 1 , have 
localization length ~ R c . Each mode can be viewed as 
harmonic oscillations of a segment of length R c in a po- 
tential well created by the pinning centers. This echoes 
the picture depicted in Fig. 2 and arguments given in 
Sec. 1. The lo 4 -\&w for the the low-frequency ("soft") 
modes follows from a careful derivation of the the sta- 
tistical distribution of the curvatures of such collective 
potential wells. One finds [22] that the density of states 
v = {8(e — e;)) of the eigenmodes goes as v(e) oc e s with 
s = 3/2. The conductivity is proportional to u, and 
once the Jacobian of the transformation from e to to is 
taken into account, one obtains Sftecr oc lo 2s+1 — to 4 . 
What is important here is that the derivation of v and 
the exponent s is rather insensitive to microscopic de- 
tails. Therefore, cj 4 -law is expected to hold in higher di- 
mensions and for other types of pinning models [24,23] . 
Shortly below, we will make use of this universality. 

Let us now return to the 2D case but first examine 
a model, which is a truncated version of the original 
one. Consider what happens if we remove by hand all 



elements of the dynamical matrix D that connect L 
and T modes. (On some crude level, it corresponds to 
a 2D crystal in zero magnetic field.) As in ID case, we 
reduce the problem to solving a Schrodinger equation. 
Now we have two such equations, 

Htut = — /iV 2 iiT + 5r(r)uT = eur (3) 

for T-modes and — W 2 ul + Sl^ul = cul for L- 
modes. Here St and Sl are the T-T and L-L com- 
ponents of the matrix S, and A should be inderstood 
as integral operator. What is the structure of the re- 
sultant eigenmodes? The T-modes with ui < uio again 
form a set of oscillators localized on the scale of the 
pinning length R c (u>o, introduced in Sec. 2, can be 
thought of as the 5 = pinning frequency). The soft 
T-modes correspond to w < coo, i.e., e <C pu>o = Ho. 
Turning to the L-modes, in principle, we also obtain a 
set of local oscillators. If A could be replaced by some q- 
independent constant much larger than p, as in the case 
of Coulomb interaction screened by a nearby metallic 
gate, then the localization length of L-modes would be 
Rl = i? c \/A//i. It is much larger than R c because Sl 
has about the same rms fluctuations as St yet the L- 
modes are much stiffer and resist localization effects 
more efficiently. If q-dispersion of A is taken into ac- 
count, then one finds [9] that the eventual localization 
does occur, but at a scale so large that it can be treated 
as infinite for all practical purposes. In this approxima- 
tion, randomness in Sl does not affect the dynamics, 
so that Sl can be replaced by So- We arrived at a sys- 
tem of localized T-phonons and delocalized L-phonons. 
In such a system, 5fte<7 has an infinitely sharp pinning 
mode. 

Let us now reinstate the L-T mixing. We will show 
that localization does occur at some reasonable scale 
Rm in this case and that Aio p is finite albeit much 
smaller than co p . However, this is not a trivial task 
because the dynamical equations become much more 
complicated. Fortunately, in the large-B limit there 
are two simplifications. First, we only need to consider 
mixing due to the Lorentz force. (The L-T mixing due 
to disorder is not important for the main results.) Sec- 
ond, the inertial terms —pco 2 in D° [Eq. (1)] can be 
neglected. The equation for ul becomes 

(-AV 2 + S )u L = ee c [H T n V. (4) 

The resolvent of the operator Ht [Eq. (3)] on the right- 
hand side of Eq. (4) can be written in terms of local- 
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ized T eigenmodes discussed above. We will see that 
Rm S> Rc', therefore, it is legitimate to coarse-grain 
the model on the scale R c , so that each localized oscil- 
lator becomes represented by a single site. This leads 
to the equation 

{-\V 2 + e )u L =^(ee c /e j )i&5(r-R j )uz,. (5) 



In the real space, it reads 



Ul(t) = ee c ^ 



K(r - Rj 



(6) 



where K(r) is the inverse Fourier transform of 
Rc/[\(q)q 2 + eo]- Function K is of the order of Kq — 
R 2 /\ at r ~ i? c . At large r, it decays as 1/r 3 . 

The spectrum of eigenmodes can be obtained from 
Eq. (6) by setting r = Rj and solving the resultant 
system of linear equations. For example, if we neglect 
randomness and replace e/s by their typical value ITo, 
we obtain e ~ IIo/e c , and so lo p ~ loq/uj c in agreement 
with Sec. 2. 

Let us now try to account for the randomness. Con- 
sider a mode with u> ~ cu p . It is easy to see that soft 
modes with energies tj <C e s = (/i/A)Ilo have a dra- 
matic effect onui. Within a distance ~ (X/fi)R c to the 
locations of such a mode, ul is strongly suppressed. We 
interpret this as an indication that these soft modes 
act as strong scatterers for the magnetophonons and 
propose that the magnetophonon localization length 
is of the order of the average distance between these 
scatterers, Rm ~ [e s ^(es)] _ (see Fig. 3). Using the 
3/2-law for the soft mode density of states we ob- 
tain Rm ~ (A//x) 5 ^ 4 _R c . Here and in previous formulas 
A = A(1/-R C ). 

Let us now evaluate Au p . To the left-hand side of 
Eq. (4) each scatterer contributes a term of the or- 
der of (A/ R 2 )(R C / Rm) 2 , where the first factor is the 
characteristic value of AV 2 in the vicinity of the scat- 
terer and the second term is the ratio of the area ~ 
R 2 where the nonuniformity in ml is concentrated to 
the total localization area R\i- The net effect of this 
term is to shift e by a random amount of the order of 
Ae ~ (X/e c )(R c / Rm) 2 ■ Such shifts produce an inho- 
mogeneously broadened pinning mode with the qual- 
ity factor Q = uj p /Auj p ~ e/Ae ~ (A//i) 3/2 > 1. This 
result was previously obtained by a somewhat differ- 
ent method in Ref. [22] where the lineshape was also 
calculated: 
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Fig. 3. Schematic structure of the localized magnetophonon 
eigenstate. Circles represent the T-modes. Soft modes that 
act as strong scatterers are shaded. 



(t(uj) = —i 



Here fi 



3/2 



rrieUJ 2 , [1 — i(ui/fl) 



3/212 _ 



(0J/UJ P ) 2 



(7) 



v/A//i. It should be mentioned that other 
theories [18,19] give different predictions for Au> p . 
Their critique is given in Ref. [9] . 



4. Comparison with experiments and 
conclusions 

The main uncertainty in comparing the outlined the- 
ory with the experiment is the nature of disorder. It 
is far from obvious, e.g., why pinning centers should 
be weak rather than strong [25] in GaAs-based 2DEG. 
However, assuming this is the case, we can take em- 
pirical values of uj p and n e and work back to the dis- 
order parameters to see if they are reasonable. In this 
manner, one obtains R c ~ 6a in samples with largest 
Q [9,7] . For the root-mean square amplitude of the ran- 
dom potential we get U ~ 0.2 K. As remarked in Sec. 1, 
the pinning frequency can also give information about 
the correlation length £ of the random potential. It is 
easy to show [18,19,9] that u p as a function of B should 
reach a maximum at the point where £ is of the order 
of the size of the bubbles (a few magnetic lengths). Us- 
ing this approach, £ of the order of a few nm can be 
estimated. Fertig [19] suggested that the disorder with 
such characteristics may originate from surface rough- 
ness. 

Regarding the linewidth, our theory predicts quality 
factors Q ~ 200 in the limit of infinite sample size 
and zero temperature. Finite width of the transmission 
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line (~ 30 /im or 10 3 a) in the experiments imposes the 
upper limit of about 30 on achievable Q. The remaining 
discrepancy may be related to quantum and thermal 
effects, which warrant further study. 
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